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^ \ Abstract 

^ , The gravitational energy-momentum within a small region as deter- 

I ^ ■ mined by two tetrad-teleparallel expressions is evaluated with the aid of 

^— i ' an orthonormal frame adapted to Riemann normal coordinates. We find 

\^ , that the gauge current "tensor" does enjoy the highly desired and rare 

£\j ■ property of being a positive multiple of the Bel-Robinson tensor, whereas 

Moller's expression does not. 
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1 Introduction: energy-momentum localization 

The localization of energy-momentum for gravitating systems is still an out- 
standing fundamental problem pQ. The classical attempts to identify a gravi- 
tational energy-momentum density for Einstein's covariant theory, general rel- 
ativity, had all led to various non-covariant expressions which could be written 
as the partial derivative of some particular non-covariant, coordinate system 
dependent superpotential, (see, e.g., [2 [3j HI [5] ) . As coordinate systems have 
! 110 physical significance, these energy-momentum density pseudotensors had no 

clear physical meaning. This led some to argue that there was no physically 
meaningful gravitational energy- momentum density, and, moreover, that this is 
■ just what we should expect from the equivalence principle (see in particular [6], 

§20.4). 

In 1961 M0ller constructed an energy-momentum density which, although 
itself still a pseudotensor, nevertheless has a superpotential which is a tensor 
under coordinate transformations [7 . M0ller achieved this "tensor" form by 
introducing an orthonormal frame, a tetrad (a.k.a. vierbein). His superpotential 
depends on the local choice of the orthonormal frame and behaves as a tensor 
with respect to coordinate transformations. Like many other energy-momentum 
expressions, the value Moller's expression assigns to a spatial region is not as 
ambiguous as one might have first thought, it is quasi-local pQ: it depends on the 
fields only at the boundary of the region. More precisely the energy-momentum 
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M0ller's tetrad expression assigns to a spacetime region depends — like other 
pseudotensors — on the boundary choice of the coordinates, but unlike the other 
pseudotensors this dependence is tensorial. Moreover, it also depends on an 
additional object which includes non-physical information, namely the choice of 
tetrad on the boundary. 

M0ller noted that his tetrad description could be given an interesting re- 
formulation in terms of teleparallel geometry. The teleparallel reformulation of 
Einstein's GR (a.k.a. the teleparallel equivalent of GR (TEGR) and GRm) has 
attracted interest not only for its presumed advantages for describing energy- 
momentum but also as a gauge theory of spacetime translations. Within the 
context of the tetrad-teleparallel theory investigators (see [51 [Ul fTUl fTTl fT^l fT5l fT^l 
n~5l [T6] and the works cited therein) have proposed another energy-momentum 
expression. It can be identified as the teleparallel translational gauge current 
density. 

Nevertheless, largely because of its perceived advantages for energy-momentum 
localization, Miller's tetrad expression (even though there is no generally ac- 
cepted frame gauge condition) — especially in its interesting teleparallel description — 
has continued to attract interest over the years (see, e.g., p!7T [18l IT9l [TOl [TTl [20] 
and the works cited therein). 

In certain special cases, however, there is a natural orthonormal frame; then 
both M0ller's expression and the gauge current yield an unambiguous energy- 
momentum. In particular this is so asymptotically — at spatial infinity. In that 
case M0ller's expression (like most others) works well (see [17] for an explicit 
verification; moreover Moller's tetrad expression in fact also works well at future 
null infinity [18j). This asymptotic success is actually not at all surprising; 
having the proper asymptotic behavior is a relatively weak requirement, for in 
this weak field region an expression need only have the proper linear theory 
limit. 

The situation is different in the one other situation where there is a natural 
frame — a case which has, to our knowledge, not been previously investigated for 
the tetrad expressions — namely the small region limit. In this limit, to zeroth 
order, one should get the material energy-momentum density — a quite weak 
requirement which follows from the equivalence principle. On the other hand 
the proposed small vacuum region limit is that, to second order, one gets a 
positive multiple of the Bel-Robinson tensor [2lJ [22j [1] (that would be sufficient 
to guarantee that the energy of a small region was positive). Now this latter 
requirement is especially interesting as a test of proposed energy-momentum 
densities, since it probes the expression beyond the linear order. It is a strong 
criterion, capable of excluding many otherwise acceptable expressions, in par- 
ticular none of the classical pseudotensors satisfy this requirement (although 
certain artificial combinations of them do [2"2l |2"5I l2"l] ). 

Here, using Riemann normal coordinates and the associated "normal" tetrad, 
we examine M0ller's expression and the gauge current in the small region limit. 
We find that the gauge current naturally satisfies this highly desirable vacuum 
Bel-Robinson property while M0ller's expression does not. 

For notation we follow 6 unless otherwise noted. Here Greek indicies are 
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used to refer to spacetime and, unless otherwise noted, a completely general 
frame. However, in those sections where it is necessary to make the distinction, 
we use Greek indicies to refer to othonormal frames, with Latin indices reserved 
for holonomic (coordinate) frames. 

2 Conserved energy-momentum densities from 
the field equations 

A gravitational energy-momentum density is easily derived from Einstein's equa- 
tions expressed in terms of differential forms: 

R a / 3 Ar ] J tl =-2KT fl . (1) 

Here k = 87rG/c 4 is the gravitational coupling constant (we will use units 
with c = 1), R a p is the curvature 2-form, T M = T v ^r] v is the source energy- 
momentum 3-form, and we are using Trautman's convenient dual form basis 
rf* " := *($ a A . . .), where is the co-frame. The left hand side of ([T]) is just 
—2G v ^it. v , the Einstein tensor expressed as a 3-form. Using the definition of the 
curvature 2-form in terms of the connection one-form and extracting an exact 
differential leads to 

R a pf\vJ» ■= (dr^ + r Q 7 a r^) a ^ 

= d(T a A rjj^) + T a A drjj^ + T Q 7 A Hp A Va % 

= d(r a p a V J\) + r* p a v\ a riaP \ - r% a n> a vJU 2 ) 

where we have used Dr] a ^ ^ = 0, which follows since the connection is metric 
compatible and torsion free. Using this expansion one can rewrite the Einstein 
equation {]]) in a neat form (which is remarkably similar to the form used by 
Einstein when he was still searching for a good gravity theory [25] ): 

dp M = 2kV^ (3) 

where the energy-momentum (superpotential) 2-form is 

:= -r% A r, a \, (4) 
and the current is the total energy-momentum density (3-form) 

P M :=t„ + T M , (5) 

which "automatically" satisfies the current conservation relation dV^ = [5]. 
This total energy-momentum current complex includes the (non-covariant) grav- 
itational energy-momentum density 

* M := (2k)- 1 (r^ a r\ a j] a \ - r% a a Va ^) . (6) 
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According to this prescription the total energy-momentum within a region is 
given by 

p*<y) ■■= I ^ = w 1 1 = (2k)- 1 I P „. (7) 

JV JV JdV 

The volume integral form would lead one to expect that the value depends 
on the quantities and choice of frame throughout the region, but the closed 
2-surface integral shows that the value is quasi-local. The value is still non- 
covariant: it depends on the choice of frame — but, as we have already pointed 
out, only on the choice at (and, through the connection, near) the boundary. 

The 2-surface integrand is 

p, := -T% A r, a \ = -T a Py ^S r ^r hp . (8) 
Expanding the components of this compact 2-form expression gives 

(r - rv/ p )<^ + (r W - r"^)^ + (r V/ 7 - r v/ T )<^- (9) 

Specializing to the case where the frame is holonomic this expression is exactly 
the superpotenial found by Freud [26 ; in that case the associated gravitational 
energy-momentum density ([6]) is the Einstein pseudotensor 3-form. On the 
other hand one can choose the frame to be orthonormal, then these same for- 
mal expressions become the those of the tetrad-teleparallel translational gauge 
current QUI EH E2 [14] , which — as we will elaborate on in the next section — 
are closely related (see [H]) to those proposed by M0ller [7] in 1961 (by the way, 
a differential form construction of these expressions virtually the same as ours 
was presented some time ago by Wallner [27], and similar arguments appear in 
|llj and [15j): the noteworthy thing is that these tetrad expression are tensors — 
under coordinate transformations. Although they are completely independent 
of the choice of coordinates (i.e., they are covariant under coordinate transfor- 
mations), they do depend on the choice of tetrad (in this important sense they 
are still non-covariant) . More specifically the energy- momentum values they 
determine are quasi-local: they depend on the choice of tetrad, but only on the 
choice at and near the boundary. 

3 M0ller's expression and the gauge current 

The traditional metric approach to gravitational energy-momentum had led 
to various pseudotensors (see, e.g., [U GH EQ E] ) , none really satisfactory. Then 
M0ller 7 replaced the metric by an orthonormal frame (a.k.a. tetrad, vierbein). 
The resulting formulation admits an interesting alternate geometric interpreta- 
tion in terms of teleparallel geometry EH EH E3 EH HI HOI E3 • It has been 
argued that this framework is more suitable for identifying a good gravitational 
energy-momentum density. Indeed, using this approach M0ller put forward his 
well-known energy-momentum expression. 

Recall that the Einstein pseudotensor can be obtained as the (Noether) 
canonical energy-momentum density from the Hilbert scalar curvature Lagrangian 
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after a certain (non-covariant) divergence has been removed (which removes all 
the second derivatives of the metric): 



£v{g,dg) := ^-S,(^ m rvC). (io) 

this is related to the aforementioned Einstein pseudotensor 3-form by t^j = 
t l E jT]j. Similarly, one can obtain M0ller's expression by using the tetrad e a i 
(related to the metric by </y = g a /3e a ie l3 j, where g = diag(— 1, +1, +1, +1) is 
the Minkowski metric) as a variable and removing an appropriate divergence 
which contains all the second derivatives of the tetrad: 

C M (e,de) := eR - d^eeje^ k 5^), (12) 

^ = = dd^k 9 ^- 5 ^- (13) 

(Here e := det e a i, the dual frame satisfies e a l e a j = 8 l - and e a l e^i = 5p, T a ^ = 
T a p{dk), and Greek and Latin indicies are transvected using respectively g a p 
and gij.) The associated M0ller 3-form is tMj = tMjVi- From this perspective 
M0ller's expression is quite natural, namely it is the (Noether) canonical energy- 
momentum density associated with the tetrad variable. It should be noted that 
exactly this same density can also be obtained from our considerations in the 
previous section — simply by formally replacing /i by j, while keeping all the 
other indices referring to the orthonormal frame. For more on these two closely 
related expressions see jTUJ [T5] . 

In sharp contrast to the metric formulation, within the tetrad/teleparallel 
formulation investigators [30[lOl[T3[Tl[T4l[2Ol[T5l[T6] have been led to only 
these two (closely related) quasi-local boundary term expressions for the energy- 
momentum within a volume V: 

P?(V) := £ Pll> T?<y) := / Pj =£ (14) 



JdV JdV JdV 

respectively, the translational gauge current and the M0ller expression f7] . M0ller 
had pointed out that his superpotential (which appears here as a 2-form inte- 
grand) is tensorial (i.e., it transforms homogeneously under a change of coordi- 
nates); however its differential, 

t M] = dpj = die^jPf,) = de^j f\p^ + e^jdp^, (15) 

the M0ller tetrad-teleparallel energy-momentum 3-form, is not a tensor with 
respect to coordinate transformations (as M0ller himself noted) — because of 
the factor de'V In contrast, it should be emphasized that both the translation 
gauge current superpotential 2-form and its differential, the gauge current 
3-form ([6]), are true tensors — under changes of coordinates. 

The tetrad theory, however, does have local Lorentz gauge freedom. The 
gauge current expressions do depend on the choice of orthonormal frame, and 
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thus still contain some observer dependent information mixed in with the phys- 
ical information in the energy-momentum expression. Nevertheless one can 
regard the gauge current expression as preferable to any of the pseudotensors 
or Moller's tetrad expression, since an orthonormal frame is more physical than 
an arbitrary choice of coordinates. 

Concerning the ambiguity re the choice of frame, it is important to note that 
the quasi-local values depend only on the choice of frame on the boundary, and 
not on the choice within the interior of the region. 

It should also be mentioned that, unfortunately, in some earlier investigations 
by our group [551 1231 US] we misidentified the gauge current as the expression 
of M0ller. (From our perspective the gauge current is the natural choice, and we 
just assumed that was what M0ller had used — without actually carefully reading 
his work. While we can appreciate that his expression is — from the Noether 
approach (|12|) . (|13|) — also a natural choice, the coordinate non-covariance of his 
energy- momentum density is certainly a liability.) 

4 Riemann normal coordinates and normal tetrads 

To find the energy-momentum within a small region surrounding a particular 
point, we look to the 3-forms V^, Vi expanding them in a power series. For 
this purpose we choose Riemann normal coordinates x l centered at the selected 
point. The Maclauren- Taylor expansion of the holonomic components of the 
metric and connection are well known (see, e.g. [BJ, §11.6): 

9ij\o = 9ij, 9k9ij\a = 0, 3d k igij\o = —Rikji — Riijk, (16) 

r 4 jfe|o = 0, 3diT z jk\o = —R l jki — R'kji- (17) 

Here g~ij — diag( — h ++) is the Minkowski metric. In the associated "normal" 
orthonormal frame, the coframe # a = e a kdx k and connection one-form T a p k alx k 
components take closely related analogous values: 

e a i\o = S a jl d k e a 3 \ o = 0, 63 fe/ e a ,| = -R a kj i - R a ijk, (18) 

r Q w | = o, 2d k r a 0l \ o = R a pkl . (19) 

It is readily verified that these values satisfy, to the appropriate order, the two 
relations which transform the metric and connection coefficients between the 
holonomic and orthonormal frames: 

9H = 9^^,, ePjVfH = T k 3l e a k - d t e a r (20) 

5 Small region values 

Here we present the energy-momentum values for small regions obtained from 
the expressions mentioned above. For non- vacuum regions all of the expressions 
reduce in zcroth order to the material energy-momentum density, in accord with 
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the equivalence principle. The value obtained for vacuum regions using the 
holonomic Einstein pseudotensor has long been known [6l [2lJ [22] . To second 
order in RNC it is 

2nVf = 2nt Ej ~ x k x l ^{AB - Sy jH riu (21) 

where 

S a 0fj,u ■— R af i\<jRf3iy Xa + Rav\oRpn Xa + -^9aP9ti^R\aK,pR X ' 7K ' P , (22) 

and 

BaPfiy ■= RaXtiaRp^ 'v a + Ra\vaRf3 X p, a — -^9at3R la& \j,R-/a8u ■ (23) 

is the celebrated Bel-Robinson tensor. (This tensor has many interesting prop- 
erties, in particular in vacuum — where the Riemannian curvature reduces to the 
Weyl curvature — it is totally symmetric and traceless.) 

5.1 The tetrad-teleparallel gauge current 

For the gauge current, expanding using Riemann normal coordinates and 
the associated normal tetrad gives, to zeroth order (unsurprisingly) only the 
source energy momentum density — just as it should according to the equivalence 
'principle. In vacuum regions "P M reduces to i M ([6]), and the leading non- vanishing 
value — using (fT9|) — appears at the second order: 



x l x m 



2«7> M = r Q/3 a r\ a Vapx - r Q 7 a a (24) 

{R aP uR X » m] - 5 x R a 7 uR^ mj ) dx l A dx j A n apx 

yri la-It f_tmS ~ O^ri ryi a JX mS) a/3XV^ 
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X' 
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4 


X' 


\ x m 



2R^ mS R vX f - ^R^^RyaSmj r\ v (25) 

. B\i mVu , (26) 

proportional to the Bel- Robinson tensor. In this calculation we have used the 
vanishing of the Ricci tensor in vacuum and some well known curvature tensor 
symmetry properties. 

To see why it is so desirable to get just the Bel-Robinson tensor one can 
integrate (|26|) over a small coordinate sphere in the surface x° — 0, using (with 
a,6,... = 1,2,3) 

x a x b d 3 x = -5 ab [ r 2 d 3 x = —S ab r 5 , (27) 
3 J 3-5 K J 

and the traceless property of the Bel- Robinson tensor to get for the gauge current 
energy-momentum 

P G ~ (2 K )- 1 B\ ab S ab ^^r 5 = B\ 00 ^r 5 = -L B \ 00 ±-r 5 . (28) 
n \ i nab 3.4.5 pUU 5! K 2G M 5! 
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This result is best appreciated when expressed in terms of the (traceless, sym- 
metric) electric and magnetic parts of the Weyl tensor, E ab := Ro a ob, Hab := 
\t a cdR cd ob- We then have a value similar to that in electrodynamics: 

Pq = (P^Pg) - £g {^{E ab E ab + H ab H ab ), e acb E ad H d b ^j ; (29) 

hence Pq satisfies an important energy condition: it is future pointing and 
non-spacelike since Pq > \Pq\ > 0. 



5.2 M0ller's expression 

Turning now to M0ller's expression Vj\ 

2k(Tj + tf) = d Pj = d(e^) = + detl j A 

= e^ i (2K)(T |t + t°) + dte^dx 1 A {-Y al3 m Ax m ) A Va ^. (30) 

To zeroth order this is again the material result one expects in accord with the 
equivalence principle. For small vacuum regions we find to lowest non- vanishing 
order 

2ntf = 2 K tf + d l e^(-r a ' 3 m )5*Xv l (31) 

T l r m 1 1 

^ — B'itaiX - + R! l n j i)x n {--R a P km x k )8%m (32) 

= ^-B i jlm r h -—(2B + S) i jlm r h (33) 

= X l x ™-L(4B-Sy ilmm . (34) 

Remarkably it turns out to be proportional to the Einstein value (|2"Tj) . (As far 
as we can see this is just an accidental coincidence.) 

According to this measure the energy within a small sphere of radius r is 

p° = ^(B--sr ab r x ^ x =^ { B- 1 -sr ab s ab ^- 

12k 4 J 12k 4 3-5 

= ^^r 5 (7E ab E ab -3H ab H ab ), (36) 
3G 5! 

which can be negative. 



6 Conclusion 

One reason that the tetrad-teleparallel formulation of GR has been favored is 
because it has been believed to have some advantage with respect to the long- 
standing problem of how to localize gravitational energy. Within this frame- 
work two energy-momentum expressions have been advocated. Here we have 
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shown that the tetrad-teleparallel gauge current (which had already been rec- 
ognized as one of the best descriptions of the gravitational energy-momentum 
for GR) satisfies another important criterion. Whereas the desired small region 
Bel-Robinson property is not satisfied by M0ller's expression, it is naturally sat- 
isfied for the tetrad-teleparallel gauge current energy-momentum density. An 
important consequence is that the gravitational energy according to the latter 
measure is positive, at least to this order. (We expected this positivity result 
since in fact there is a positivity proof for the energy associated with the tetrad 
gauge current expression [5].) 

We stress that the vacuum small region Bel-Robinson property is, as ex- 
emplified by the two cases considered here, a strong test capable of excluding 
many otherwise acceptable expressions; indeed none of the classical pseudoten- 
sors (and in this category one can include M0ller's tetrad-teleparallel expression) 
satisfies this requirement (although certain quite artificial combinations of them 

do [221 EH EH). 

Compared to these the tetrad-teleparallel gauge current energy-momentum 
density stands out. It is certainly a better description for gravitational energy- 
momentum. In addition to being a tensor (under coordinate transformations) it 
also enjoys the highly desired and rare property of having its small region value 
be a positive multiple of the Bel- Robinson tensor. 
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